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PRELIMINARY SURVEY CF PROPULSION USING CHEMICAL 
ENERGY STORED IN THE UPPER ATMOSPHERE 
By Lionel V. Baldwin and Perry L. Blackshear 


SUMMARY 

This report presents a preliminary study of a ramjet that would use 
the chemical energy of dissociated molecules in the ionosphere for pro- 
■ H puls ion . A review of the physical properties and chemical composition 

^ of the upper atmosphere shows that the available energy is not sufficient 

o for flight requiring aerodynamic lift. Above 300,000 feet, dissociation 

energy might be useful for satellite -sustaining. Comparison of maximum 
thrust and external drag for an orbiting ramjet from 300,000 to 700,000 
feet shows that a positive net-thrust margin might be obtained over the 
entire region for seme of the engine sizes considered. However, the re- 
gion between 325,000 and 400,000 feet is the most favorable, because 
high gas densities produce high recombination rates. 

Two thermodynamic engine cycles are considered for a ramjet orbiting 
at 328,080 feet: (l) all- supersonic flow with frozen composition in the 

diffuser and chemical equilibrium in the nozzle j (2) normal- shock inlet 
with equilibrium expansion. Only the all-supersonic cycle shows promise 
for sustaining a satellite. Calculations of the rate of recombination 
are only approx ima te because of uncertainties in the data on the chemical 
reaction rates. These approximate calculations indicate that, although 
equilibrium expansion cannot be obtained, up to 50-percent recombination 
may be possible for a nozzle 100 feet long. This would give thrust in 
excess of the drag. 

Aerodynamic -heating problems of the orbiting ramjet axe considered. 
External surfaces are sufficiently cooled by radiation. The internal 
throat area, which presents the most difficult cooling problem, requires 
cooling rates that are met by current technology; but, in order to mini- 
mize heat loss due to internal cooling, the inlet radius of the engine 
will probably have to be greater than 20 feet for an engine 100 feet 
long. 


The recombination ramjet appears to be a marginal device even in 
the optimistic view. Before definite conclusions can be drawn, the 
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atmospheric condos it ion near 300,000 feet and the recombination-rate 
coefficient for atonic oxygen must be accurately known. If favorable 
composition and rate coefficients are measured in future research, then 
the "fuelless" recombination ramjet can probably be developed, but only 
for use where most of the lift comes from the orbiting velocity. 


introduction 

About 75 years ago, W. N. Hartley observed the abrupt termination 
of the ultraviolet end of the spectrum of all heavenly bodies and postu- 
lated the existence of ozone in the Earth's atmosphere to account for 
this absorption. In 1930, Sydney Chapman showed theoretically that the 
absorption of solar ultraviolet light should^cause diatomic oxygen to 
dissociate at altitudes above 100 kilometers. The theory was clear-cut 
and was universally accepted, although the atomic oxygen in the iono- 
sphere had not been experimentally measured. More recent theories have 
postulated that nitrogen is also dissociated in the upper regions of the 
ionosphere. Balloon experiments and, in the last ten years, high- 
altitude rockets have proved these speculations to be correct. 

A small fraction of the solar energy absorbed by the Earth is locked 
in reactive chemical species such as ozone, dissociated molecules, and 
gaseous ions in the upper atmosphere. This energy may be available for 
flight propuls ion in the upper atmosphere. For example, in 1950, Proell 
(ref. l) proposed the use of "fuelless box-ramjets" using the energy of 
atomic oxygen at about 100 kilometers to move and direct space stations, 
alter orbits, and overcome frictional losses. The "fuelless" feature of 
an engine that converts the chemical energy of reactive air particles 
into thrust is very attractive for these applications. Furthermore, the 
simplicity of the ramjet is desirable where long service times are 
required . 

The availability of atmospheric energy for propulsion is investigated 
in detail in this report in the light of atmospheric data obtained by 
rockets in the last five years. To facilitate the ramjet cycle calcula- 
tions, the reaction enthalpy for the conversion of photochemically dis- 
sociated air to stable molecules in equilibrium is summarized as an 
available energy density (Btu/cu ft of air) for altitudes from 50,000 to 
1,000,000 feet. An analysis comparing the available energy with the 
energy required for flight shows that sufficient thrust can be provided 
only by an orbiting engine . 

Next, the orbiting ionosphere ramjet is considered in detail. The 
external drag for various nacelle gecmetries_ is compared with the maximum 
available thrust from 325,000 to 700,000 feet. One of the more premising 
flight altitudes and configurations is chosen for a thermodynamic cycle 
calculation. Finally, the problems of recombination chemical kinetics 
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and internal and external heat transfer are examined for this ramjet 
orbiting at 328,080 feet. The calculation of heat transfer and friction 
drag is discussed in appendix D by James F. S chmi dt - 

Unfortunately, existing data are not complete enough to evaluate 
definitely the recombination ramjet. The purpose of this report is to 
present a prel iminar y study and to point out the research data that are 
necessary. By making plausible assumptions where data are lacking, it 
is possible to propose the major design requirements and to indicate the 
probable magnitude of the problems likely to be encountered. 


MODELS OF UPPER ATMOSPHERE 
Physical Properties 

The pressure, density, temperature, and composition of the atmos- 
phere extrapolated up to an altitude of 1,600,000 feet are s ummar ized 
in engineering units in tables I to III. Three data summaries of upper- 
atmospheric rocket research (refs. 2, 3, and 4) have been selected to 
serve as "standard atmospheres" for this discussion. The combined re- 
sults of all rocket flights from. 1946 to January 1952, presented by the 
Rocket Panel (ref. 2), are summarized in table I. Another summary (ref. 
3) was prepared by Dr. E. 0. Hulburt, Director of Research at the Naval 
Research Laboratory. Hulburt presented smoothed, average values obtained 
by a ll ,, flights up to July 1954, and he extrapolated the data from. 

725,000 to 1,600,000 feet. This working model of the atmosphere is pre- 
sented in table II. A more recent standard atmosphere (1956) has been 
proposed as the Air Force ARDC Model Atmosphere (ref. 4). Table HI is 
an abridged form of this model. 

Symbols are defined in appendix A, and the reliability of the tabu- 
lated physical properties of the atmospheric models is discussed briefly 
in appendix B. 


Chemical Compos it ion 

The chemical compositions in tables I to III are less certain than 
the physical properties, because composition data above 250, 000 feet are 
very scant. In fact, the compositions in table I are hypothetical; the 
Rocket Panel assumed these values only to calculate an average molecular 
weight, which in turn could be used to calculate temperature from, meas- 
ured pressure and density. As discussed in appendix B, the Hulburt-NRL 
model compositions from 357,000 to 420,000 feet (table H) are based on 
data from three Aerobee rocket flights in 1953. Above 426,000 feet, the 
compositions are based on assumed atmospheric conditions. The ARDC com- 
positions of table ITT were prepared from smoothed individual particle 
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concent rat ions read from a graph published in reference 5. Therefore, 
table III compositions are approximate, but they probably represent the 
ARDC adopted values within 3 percent for any component . Appendix B in- 
cludes a discussion of the compositions of the three model atmospheres. 

The chemical energy associated with the nonequilibrium composition 
of "air" at these high altitudes has been estimated frcm. the tabulated 
composition data of tables I and II (fig. l). The available chemical 
energy is taken to be the reaction enthalpy for the conversion of the 
photochemically dissociated air to stable molecules in chemical equilib- 
rium. Chemical propulsion schemes based on the presence of ozone, atomic 
oxygen, and atomic nitrogen can be evaluated readily from figure 1, which 
shows "available chemical-energy density" (Btu/cu ft of air) from 50,000 
to 1,000,000 feet. Table III energy density, though not shown, would be 
between the Rocket Panel and Hulburt-NRL curves in the region from about 

300.000 to 600,000 feet. 

The calculation procedure used to obtain figure 1 from tables I and 
II, together with seme additional composition data on ozone (table IV), 
is included in appendix B. Therefore, only a qualitative description of 
figure 1 will be given here. The solid curve between 50,000 and 260,000 
feet represents the chemical energy of ozone, which is not shown in 
either table I or II because it is present in only minute concentrations. 
The vertical line at 270,000 feet emphasizes the uncertainty of the Rocket 
Panel's assumed composition at this altitudej actually, the abrupt rise 
in atomic oxygen concentration listed in table I would lead to a nearly 
vertical slope. The solid curve frcm 270,000 to 394,000 feet is the en- 
ergy d ens ity associated with atomic oxygen; above 394,000 feet, atomic 
nitrogen begins to contribute with increasing importance. 

For the energy-density curve from the Hulburt-NRL model of table II, 
the atomic oxygen begins to contribute to the energy density at about 

150.000 feet. The dashed curve from 150,000 to 325,000 feet reflects a 
rapid buildup of atomic oxygen concentration which, though below 2 per- 
cent of the total air, is still signif icant_in terms of energy because 

of the relatively high total air density. Because the atomic oxygen con- 
centrations are small in this transition region from 150,000 to 325,000 
feet, the percentages are not listed in table II but are included in 
table TV instead. The peak in the dashed curve (fig. l) at about 

320.000 feet is caused by the leveling off of the initially steep in- 
crease in percentage of atomic oxygen in the transition region to the 
values shown in table II above 328,000 feet. The slope of this curve 
following the peak reflects the exponential decrease in total air den- _ 
sity. Atomic nitrogen starts to contribute to the energy at about 

460.000 feet, causing the change in slope shown in figure 1. 
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THE IONOSPHERE RAMJET 

Assuming that the atmospheric energy described in the previous sec- 
tion is available for propulsion, it is important next to deter min e the 
flight speeds and altitudes at -which this energy would be sufficient. 

Such an orienting analysis is presented in appendix C, which formalizes 
the relations concerning required energy, lift-drag ratios, engine effi- 
ciency, friction drag, altitude, and flight speed. This analysis shows 
that for a supersonic ramjet the available energy would not be useful 
for sustaining flight requiring aerodynamic lift. However, for an orbit- 
ing ramjet the available energy may be sufficient to overcome friction 
drag above about 300,000 feet. 

Therefore, an orbiting ionosphere ramjet is considered in this sec- 
tion. Appearance, operation, and design problems are discussed. Ho 
attempt is made to consider all possible engine designs . Rather, one 
plausible ramjet is studied in detail to test the feasibility of using 
atmospheric energy. The external friction drag for various truncated- 
cone nacelle geometries and orbit altitudes is presented first. Then, 
far a promising configuration and altitude, a thermodynamic cycle is 
followed through to determine some factors Influencing engine efficiency. 
Since the chemical kinetics for the recombination of dissociated air 
will govern the actual engine cycle, the assumptions used in the thermo- 
•* dynamic analysis are examined with respect to the kinetics. Finally, 

since aerodynamic heating can destroy an object orbiting in the atmos- 
phere, the internal and external heat-transfer problems are considered. 


Thrust-Drag Comparison at Various Altitudes 

The thrust F of a recombination ramjet is easily related to the 
energy density e of the air at the orbit altitude. The thrust per 
unit inlet area Aj_ is 

j- « nei (1) 

A 1 

The thrust is a function of engine efficiency T| and altitude, since 
e depends on altitude. If the engine efficiency is assumed to be 100 
percent, then it is convenient to define a thrust parameter that is a 
function only of orbit altitude: 

F = eJ 

A lPl Pi 

The thrust parameter has been evaluated for altitudes frcaa 300,000 to 
700,000 feet using the two model atmospheres of figure 1. The two 


( 2 ) 
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curves of figure 2 represent the probable upper and lower limits on the 
energy available for propulsion for a recombination ramjet at orbital v 

velocity . 

A similar term can be defined for external drag, but the nacelle 
geometry must also be specified. Figure 3 shows the simple truncated- 
cone nacelle geometries considered in this preliminary study. Nacelles 
having these shapes can be distinguished for the following discussion 
by giving the ratio of length to capture radius the cone half- 

angle 9, and the over-all length A variety of combinations of the 
three variables was chosen, and the external drags for each nacelle were 
calculated over an altitude range. In the friction-drag calculations 
outlined in appendix D, a procedure described by Eckert (refs. 6 and 7) 
was used. The external drag is divided by the inlet area and ambient 
density for direct comparison with the maximum thrust. 

The energy available is compared with the energy required to over- 
come drag for altitudes from 300,000 to 700,000 feet and for Jj? of 100 
and 1000 feet in figures 4(a) and (b), respectively. Figure 4(b) is in- 
tended only to show the effect of engine size an external drag and max- 
imum thrust over a wide range . Both plots give the drag parameter for 
S ?/ % of 2.5, 5, 10, and 30, and for 9 of 2°, 0°, -4.3°, and -8.6°. 

The requirement for sustained satellite flight is i> 

_JL_ > ( 3 ) 

A lPl ' A lPl W 

Figure 4 indicates that, even with an engine efficiency of 100 percent, 
the smallest engine size shown (jgf = 100 feet) would not sustain satel- 
lite flight unless very small and negative 9 could be used. As 

might be expected, larger inlets give increasingly favorable thrust mar- 
gins with this idealized comparison. 

The difference between the thrust and drag parameters in figure 4 
is proportional to the maximum net thrust. Between altitudes of 320,000 
and 400,000 feet, four of the design configurations show positive net- 
thrust capabilities. It will be shown later that the chemical recombi- 
nation reaction is strongly influenced by air density, being more favor- 
able as density increases. Furthermore, the drags in figure 4 are prob- 
ably too low above 450,000 feet, because the flow would be approaching 
free -molecule flow (see appendix C). Therefore, for the following cycle 
analysis, 328,080 feet was chosen as a favorable altitude, although this 
precise choice was arbitrary. 
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Thermodynamic Cycle for Recombination Ramjet 

The engine efficiency of a recombination ramjet will certainly be 
less than 100 percent , and the nacelle geometry can be expected to affect 
efficiency as the ratio of inlet to exhaust area (and therefore the de- 
gree of expansion) is varied. These performance variables an d others 
are discussed in the following paragraphs . 

The problem is to get the energetic ambient air on board, to recom- 
bine the atomic oxygen, to exhaust the hot diatomic air, and thus to ob- 
tain thrust. A circular orbit of 328,080 feet and the atmosphere of the 
Rocket Panel (table X) are the assumed flight conditions . Two different 
cycles are analyzed. However, it should not be inferred that these in- 
clude all possible cases. 

All- supersonic cycle . - The first cycle is an all- supersonic cycle; 
that is, the flow throughout the engine is above sonic velocity. The 
cycle can be divided into three stages, as follows: (l) frozen- 

composition compression, (2) conversion of compressed air to chemical 
equilibrium, and (3) chemical-equilibrium exhaust expansion to ambient 
pressure. There is an apparent paradox in these assumptions, because 
the first requires the recombination rate to be zero, while the third 
assumes that the rate is fast compared with flow residence times. For 
the moment, assume that each limiting condition could be approached by 
making the diffuser residence time short compared with the exhaust- 
nozzle flow time. Later, a kinetic analysis will estimate the accuracy 
of these limiting assumptions. 

Unfortunately, the published thermodynamic properties for air do not 
extend over the entire low-pressure region of this cycle (e.g., ref. 8). 
However, a report on the thermodynamic properties of the three -component 
hydrogen-oxygen-nitrogen system has recently been published (ref. 9). 

This system is available at the HACA Lewis laboratory in a digital - 
ccmputor program designed to study rocket nozzle perfor man ce - The simple 
expedient of driving the hydrogen content to negligible proportions 
(lO - -^ percent or less) converted the hydrogen- oxygen-nitrogen system of 
reference 9 to an argon-free air system that included the entire pres- 
sure range of interest here. All thermodynamic properties used in this 
section were obtained from this program and reference 9, using "air" 
with 1 mole of oxygen (32.000) for every 3.773 moles of nitrogen 
(28.016). 

Figure 5(a) shows the thermodynamic cycle on a plot of static tem- 
perature against static pressure. Figure 5(b) gives the corresponding 
air enthalpy as a function cf static pressure. The three stages of this 
cycle axe discussed briefly in the following paragraphs . 



8 


MCA IN 4267 


For the 328,080-foot orbit, the frozen composition is 19 percent 0 
(by volume), 9.5 percent 0 2 , and 71.5 percent Ng. The enthalpy at 

any temperature station in the c empress ion is given by 

% = [o.19(^ q ) t + 0.095(jT 02 ) t + 0.715^)^ ^TH ( 4 > 

The air enthalpy of figure 5(b) was calculated from equation (4) and the 
molar enthalpies tabulated in reference 9 for all T. The rela- 

tion between temperature and pressure for the frozen compression was 
assumed to be that for isentropic flow: 



Since the exponent r/v-l is a function cf temperature, a small temper- 
ature difference (Tjj " “ 180° R) was used, and the pressure ratio was 

calculated for these temperature increments. Reference 9 data on Cp/^P 
for the individual chemical components were averaged on a molar basis 
(as in eq. (4)) for use as y/y-1. Therefore, the compression in both 
enthalpy-pressure and temperature -pres sure relations includes an allow- 
ance for adjustment in specific heats. 

The problem of where to terminate the compression is unique to the 
recombination-ramjet cycle. The procedure used to calculate the 
compression-termination station, where the conversion of frozen to equi- 
librium is assumed to occur, is illustrated by figure 5(c). This figure 
is simply an enlarged temperature-pressure diagram, in the region near 
station 2 of figure 5(a) . 

If the frozen compression is assumed to persist until very high 
static temperature is attained, then the conversion of the ambient com- 
position air to equilibrium will lead to more dissociation and a temper- 
ature drop. To illu strate, consider station b (fig. 5(c)), which is an 
arbitrary station in the frozen compression. The constant-pressure 
adiabatic conversion from frozen to equilibrium causes a temperature 
drop to T-jj. On the other hand, let the compression be terminated at a 

low pressure (station a) . Then recombination occurs when the gas is 
converted at constant pressure to equilibrium, and a temperature rise to 
T a is noted. The compress ion-terminat ion station 2 in figures 5(a) and 

(b) is the adiabatic, isobaric, isothermal conversion of the frozen in- 
let air to equilibrium. 

It is not apparent from the preceding discussion that station 2 is 
the most advantageous conversion point, and indeed it might not be. 
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However, of the three points discussed, the isothermal conversion was 
observed to give higher engine efficiency with full expansion than either 
station a or station b. This efficiency difference was small, being less 
than 1 percent for all three cases. Therefore, the choice of any other 
conversion station in the region shown in figure 5(c) rather than sta- 
tion 2 would have had negligible effect on the results of the cycle 
analysis . 

The final stage of the cycle is an isentrcpic expansion in chemical 
equilibrium from, station 2 to station 3. In figure 5(a), note the unu- 
sual temperature-pressure dependence from, station 2 down to about 10“ p 
pound per square foot; in this region the recombination of atomic oxygen 
is occurring. The release of the chemical energy is also shown clearly 
on the enthalpy-pressure plot (fig. 5(b)). 

The chemical composition during the equilibrium expansion from sta- 
tions 2 to 3 is given in figure 6. At station 2, the difference between 
the inl et frozen composition and the composition of the equilibrium mix- 
ture is due solely to 2.23 percent HO. The equilibrium compositions of 
figure 6 are used in the chemical kinetics analysis of the next section. 

The ideal engine efficiency and the net thrust were calculated from 
the cycle summarized in figures 5 and 6. The flow conditions at the 
throat (station 2) are as follows: 

0 2 = A R~- - S 1 )2g (; J (6) 


where 


U-l = 26,050 ft/sec (orbital velocity) 
= 984.8 Btu/lb(M) 1 
Hg = 2309.4 Btu/lb(M)J 
U 2 = 24,743 ft/sec 


fig. 5(b) 


u 2 

Mo = - ■ = 7.5 

Vr^gcT 2 


( 7 ) 


and the ratio of inlet to throat area is calculated from continuity: 

Aq = U2^1^2^1 

h U 1 T 2 P 1 J( 2 


= 740.6 
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Ihe average molecular weights J£ at stations 1 and 2 were calculated 
from compos it ions given in table I and figure 6. 


In a similar manner, the ratios of throat to exit velocity and throat 
to exit area were evaluated at various stations throughout the equilibrium 
expansion. Then, the thrust per unit inlet area was determined from 


JL _ p i u i 

A 1 " 8 C 


C U 3 - U l) + P 3 




( 9 ) 


That is, the thrust depends on the inlet area and the ratio of inlet to 
exit area. If truncated-cone nacelles are used, then the drag also de- 
pends on the ratio of inlet to exit area. This suggests that an optimum 
area ratio exists, as will be discussed shortly. 

Closely related to the thrust is the engine efficiency: 


' ’ Cio) 

The Rocket Panel energy density e from figure 1 at 328,080 feet was 
substituted into equation (lO) together with the results from equation 
(9). The resulting engine efficiency is given in figure 7 as a function 
Of A-j/Aj. 

If the thermodynamic cycle assumptions of equilibrium expansion and 
reversible adiabatic flow are good approximations, then an optimum con- 
figuration for maximum net thrust can be found. Figure 8 gives the 
thrust and drag as a function of the ratio of inlet to exit area for a 
100-foot engine having an inlet radius of 20 feet. The drag is much 
more sensitive to this area ratio than the ideal thrust, and it is clear 
that the difference between thrust and drag will pass through a maximum. 
Figure 9 shows the net thrust for engines 100 feet long with various in- 
let radii as a function of ratio of inlet to exit area. The large effect 
of inlet radius could be expected, because the thrust is directly propor- 

p 

tional to the inlet area, or (R-^) , while the drag is proportional to the 
nacelle surface area, or approximately Rj_. The optimum net thrust peaks 
more sharply for large inlets, but in all cases shown the peak occurs at 
area ratios between 2 and 6j the corresponding cane half-angles are neg- 
ative, about -4° to -6°. 

Normal-shock- inlet cycle . - The cycle discussed has assumed that 
the flow in the orbiting ramjet is supersonic throughout. Many alter- 
native cycles involving shock inlets are possible, of course, but only 
the extreme case of a normal- shock- inlet cycle has been examined for 
this preliminary study. Appendix E gives an analysis of a recombination 
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ramjet -with, a normal- shock inlet under the identical flight conditions 
i of the supersonic analysis. This analysis shows that the extremely high 

static temperatures after shock result in more dissociation t han can "be 
regained in a full equilibrium expansion. A ramjet with a normal- shock 
inlet in a 328,080-foot orbit failed to produce useful thrust by at 
least an order of magnitude. 


Chemical Kinetics of Recombination 
o 

t- The recombination ramjet could supply sufficient thrust to sustain 

itself indefinitely in circular ionosphere orbits under ideal circum- 
stances: reversible, adiabatic flow and chemical-equilibrium nozzle ex- 

pansion. In this section the assumptions of frozen- compos it ion ccmpres- 
sion and chemical-equilibrium, nozzle flow are examined; the next section 
is concerned with the heat-transfer problem. 

If the chemical rate processes were understood completely and all 
kinetic data were available, the problem of adiabatic nozzle expansion 
with chemical reactions would be solvable, at least In principle. How- 
ever, even for the simple case of recombining atomic oxygen in the ab- 
sence of foreign fuels, plausible assumptions must replace data that are 
lacking today. The recombination kinetics are analyzed here for two 
limiting conditions to place upper and lower limits on the probable de- 
gree of recombination . The limiting cases considered are: 

(1) Near-equilibrium flow: This condition would occur if the re- 

combination of the atomic oxygen were sufficiently rapid to keep pace 
with the shift of equilibrium compositions resulting from, the drop in 
static temperature and pressure in the nozzle flow. 

(2) Near-frozen-composition flow: In this case the recombination 

kinetics are so slow that negligible recombination occurs . 

The details of the kinetic analysis are deferred to appendix F, 
where the general method described in reference 10 is applied to re- 
combination of atomic oxygen in the exhaust nozzle. However, a quali- 
tative understanding of the analysis is sufficient to make the results 
meaningful. In an actual nozzle, the flow may either be in chemical 
equilibrium, have frozen composition, or vary anywhere between these 
limiting cases depending on the reaction rate. For example, consider 
a hypothetical nozzle with a very large pressure ratio. The nozzle flow 
could follow chemical equilibrium initially. As the recombination pro- 
ceeds, the concentrations of atomic oxygen reactants would drop to such 
low values that the reaction rates could no longer keep pace with the 
m-tn-tmTrm rate required for near-equilibrium flow. The reaction would 
proceed in the nozzle with rapidly decreasing rate as the atomic oxygen 
was converted to molecular oxygen; but eventually the rate would become 
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60 slow that, for all practical purposes, no__more recombination would 
occur. At the latter condition, the flow is at the other limit, near- 
frozen-camposition flow. The analysis of the nozzle flow given in appen- 
dix F makes possible the numerical evaluation of the minimum, concentra- 
tion of atonic oxygen that would exist in a given nozzle when the flew 
deviates from equilibrium because of insufficient reaction rate. The 
analysis also places an upper limit on the atomic oxygen c cncentration 
corresponding to a reaction rate which, if exceeded, would cause the 
flow to deviate from near- frozen composition. 

The numerical example summarized in table V is for the nozzle ex- 
pansion discussed in the supersonic cycle. The first c olumn lists 

recombination-rate coefficients kp from 10 11 to lO^ (g-mole/liter)“ 2 
sec“-*-j cgs units are used here for convenience, because these are the 
usual units used by physical chemists reporting kinetic data. Only one 
value of kj> would be listed if rate data were available ; hut, as will 

be discussed shortly, there is enough uncertainty concerning kg to 
warrant the wide range of magnitudes in the table. 

The second column in table V gives the minimum concentration of 
atomic oxygen ( ^q, g-mole o/liter) required for near- equilibrium flow 

with the corresponding rate coefficient. Atomic oxygen concentrations 
smaller than those shown would cause the flow to deviate toward frozen 
flow. However, atomic oxygen concentration depends not only on the mass 
fraction of atomic oxygen but also on the air density, and both vary in 
the nozzle flow. That is, the minimum concentration necessary for near- 
equilibrium flow is meaningless as a guide to recombination, unless it 
is related to the temperature, pressure, and mean molecular weight of 
the expanding air. If it is recalled that the air is assumed to follow 
chemical equilibrium until the minimum concentration for near-equilibrium 
is reached, then it follows that the equilibrium concentrations in fig- 
ure 6 together with the temperature -pres sure relation in figure 5(a) give 
the required relation between temperature, pressure, mean molecular 
weight, and the mass fraction of atomic oxygen. For equilibrium flow, 
appendix F shows that the mass percentage of atomic oxygen recombination 
can be calculated as a function of the atomic oxygen concentration in 
the nozzle flow. Figure 10 is a plot of this relation. The minimum 

values in table V are related by figure 10 to the mass percentage 

of the atomic oxygen that has recombined to molecular oxygen; this num- 
ber is significant, since it sets a lower limit on the amount of recom- 
bination attained if the rate coefficient is kg. 

The last two columns of table V give similar results for near- 
frozen flow in the nozzle. If 9* 0 is smaller than the listed values 

at the given kg., then negligible reccanbinatian will occur. The relatiah 
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between <£q and mass percentage recombined cannot be determined rigor- 
ously for the near-frozen process by this simple analysis, because the 
necessary flow relations are unknown once the flow varies from equilib- 
rium. However, by assuming the equilibrium flow relation between 

and mass percentage recombined (fig. 10} , it is possible to set an upper 
limit on the degree of recombination. Regardless of what the actual re- 
action rate is in the intermediate flow between equilibrium and frozen, 
it will certainly be slower than in equilibrium flow. Therefore, the 
corresponding mass percentage recombined given in the table was taken 
from figure 10, and it is the maximum degree of completion possible in 
the nozzle for a particular value of k^. 

A recombination-rate coefficient of 10"^ (g-mole/literj^sec - ^ is 
frequently employed by physical chemists for reactions of this type, be- 
cause, although no experimental data for atomic oxygen are available, 

values of k^ of 10 10 have been observed for iodine, bromine, deuterium, 

and hydrogen atomic recombinations (ref. ll). Furthermore, the reaction 
rates were found to be almost temperature- independent, with activation 
energies near zero. This is the only justification for referring to 
single-valued k^ in a nozzle where the temperature varies over a 4000° 

F range. If lO 1 - were correct, then the flow would never approach chem- 
ical equilibrium, and 8-percent recombination would occur at most. Under 
these circumstances, the ionosphere ramjet would not be feasible. 

Fortunately, there is reason to believe that the rate constant for 
recombination of atomic oxygen may be much larger than 10^-. Feldman 
(ref. 12) has recently reported an interesting high-temperature shock- 
tube study of the kinetics of "air." Air was dissociated at the high 
temperatures following the passage of the shock through the tube, and a 
preliminary experiment established that the dissociated air whb in chem- 
ical equilibrium. The hot equilibrium air was then passed over a re- 
versed wedge in a Prsndtl-Meyer expansion. By observing Mach angles in 
the expansion, Feldman was able to conclude that the air bad maintained 
chemical equilibrium throughout the rapid drop in static temperature and 
pressure . By using an analysis similar to that given in appendix F, he 
calculated lower-limit recombination-rate coefficients kp from 4500° 

o v 

to 7900 K and from 0.005 to 1 atmosphere for the following reaction: 

kf 

s/ + + X ^ + X 

^r 

where represents a molecule of "air" and an atcm of "air" 

having half the molecular weight of the molecule. Lower-limit coeffi- 
cients kp as high as 10"^ (g-mole /liter ) _ ^sec _ - were reported. 
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If 10 13 is the correct kp, then table V Indicates that, although 
near-equilibrium flow would not be attained anywhere in the nozzle, over 
50-percent recombination is possible in a nozzle 100 feet long. The re- 
sults of the AVCO shock- tube research can probably be applied to t he re- 
combination of oxygen, because atomic oxygen is the predominant dissoci- 
ation species of "air" in the temperature -pressure region where 10^ was 
deduced . 

The effect of changing the length of the adiabatic nozzle can also 
be obtained from table V. The values shown are for a 100-foot nozzle, 
hut it is reasonable to expect that going to longer nozzles will give 
more time for recombination. Equations (F16) and (F2l) of appendix F 
give this quantitative statement. For example, if the nozzle were 1000 
feet long, then the effect would be to divide each kf in table V by 10. 
Thus, for k f = 10 13 , up to 78-percent completion might be expected, but 
equilibrium flow would not be approached anywhere in the nozzle. 

The preceding discussion of table V assumes that the inlet- diffuser 
flow will have frozen composition. If recombination -rate coefficients 
of 10l 3 or less are correct, then the assumption of frozen inlet flow is 
probably much better than the assumption of equilibrium nozzle flow. 

This is especially true where the inlet diffuser length is about 10 feet 
compared with at least 100 feet for the exhaust nozzle. The major prob- 
lem will be obtaining sufficient recombination . Any recombination occur- 
ring in a short diffuser will simply reduce both the design temperature 
and the pressure at station 2. That is, the limiting cycle of figure 
5(a) probably envelops the actual cycle, which would have (l) a slightly 
steeper temperature-pressure dependence toward the end of the compres- 
sion, (2) a lower tenperature -pressure design station for the start of 
the nozzle expansion, and (3) a steeper temperature-pressure curve in 
expansion with less pronouncedly flat recombination region. Katurally, 
the enthalpy difference between stations 1 and 3 would be less than the 
limiting cycle, because, unlike that in figure 5, the exhaust gas would 
probably have some dissociation energy. The actual thrust would be less 
than that calculated in this cycle analysis. 

A kinetic solution to give values for the actual cycle outlined in 
general will, not be attempted here. The uncertainties in both atmos- 
pheric composition and reaction-rate coefficient are much too large to 
justify laborious "exact" numerical designs. However, if about 50- 
percent conversion could be obtained in an engine about 100 feet long, 
then by making the capture radius Rq_ sufficiently large, the recom- 
bination ramjet principle would be worthy of further study. 

The chemical kinetics of recombination are very sensitive to orbit 
altitude, since the reaction rate is roughly proportional to the cube of 
the ambient density. The results in table V are applicable to an 
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altitude of 328,080 feet in the atmosphere of table I, which is the low- 
est region of appreciable photochemical dissociation in this model. 

Lower orbits, though desirable kinetically, appear to be Imp ossible. 


Heat-Transfer and Cooling Problems of Orbiting Ramjet 

External . - The external equilibrium wall temperature for the re- 
verse truncated- cone nacelle is radiation-controlled. Figure 11 gives 
the results of a calculation of the wall temperature at the 100-foot 
point on a -8.6° nacelle over an altitude range from. 328,000 to 750,000 
feet. The calculation assumed that the surface of the nacelle was normal 
to the sun's rays and that it radiated energy as a black body to space 
(see appendix D) . Within the accuracy of the calculation, the convec- 
tive heat transfer did not contribute at = 100 feet; but, for cone 
half -angles more positive than -8.6°, a small variation in the wall tem- 
perature was noted over this altitude range. The variation of wall tem- 
perature over the length of the nacelle is given by figure 12. An orbit 
altitude of 328,080 feet was chosen for this example, and the wall tem- 
peratures for cone half -angles of -4.3° and -8.6° are typical of the 
effect of convective heat transfer on the surface linear dimension. The 
important conclusion to be drawn from figures 11 and 12 is that, because 
of radiation cooling, external heating will not be a problem. 

Internal . - The internal surface, having no radiant cooling outlet, 
presents the major cooling problem of the orbiting ramjet. The throat 
of the engine (station 2) is most difficult to cool, because the maximum 
heat-transfer rates occur there. To estimate the magnitude of the prob- 
lem, consider a ramjet with R-j_ of 20 feet orbiting at an altitude of 

328,080 feet. The ratio of inlet to throat area A-^/Ag frcm. the cycle 
calculations is 740.6. 


The heat-transfer coefficient is calculated frcm station 2 condi- 
tions and the Dittus-Boelter equation, which is commonly used in rocket 
throat cooling problems (ref . 13 ) . As shown in appendix D, the maximum 
heat -transfer rate for the throat of this ramjet is 



10.82 


Btu 

(sec)(sq ft) 


Current rocket practice for throat cooling can handle heating rates as 
great as 100 times 10.82. Therefore, the hot-spot cooling for the or- 
biting recombination ramjet could be solved with present technology, 
the problem being less difficult than that encountered in rocket design. 

Of course, the entire internal surface of the engine will require 
cooling. The cooling system might use the outer nacelle as a radiator 



16 


NACA TN 4267 


and use a coolant in a continuous cycle. The heat transferred "by the 
radiator could make the ramjet inoperative if this heat were a large 
fraction of recombination energy, because this radiated energy is not 
available for propulsion. The fraction of the heat lost in cooling will 
depend on the geometry of the engine and the orbit altitude, but it is 
interesting to make a crude analysis at 328,080 feet to complete this 
preliminary survey. 

The heat-transfer rate at any section of the internal surface is 

■TjjjL- = h^(T aw - T w ) (ll) 


The adiabatic wall temperature will vary only slightly in the engine at 
hypersonic speeds > therefore, if it is assumed that the wall temperature 
is maintained constant throughout, T &w - Ty. does not vary in the engine. 

By substituting the Dittus-Boelter equation for the heat-transfer coeffi- 
cient (appendix D, eq. (D19)), equation (ll) becomes 



The velocity is nearly constant in hypersonic flow; therefore, from con- 
tinuity the following approximations can be made: 

U «* U = 26,000 ft/sec (l3) 

pd 2 - p-jd 2 (14) 


Equation (12) can be written in terms of equations (13) and (l4) for the 
total cooling required: 


^cooling 


/ 

-\0.8 

0.023 *yPr^ ,33 (- 

— ) (T 
LjcJ K aw 


-T> 


0.8 

1 




(15) 


The numerical value of equation (15) depends strongly on the inlet di- 
ameter d-j_ and to a lesser extent on the internal design (i.e., 

d = f (x) ) . The greatest value of q 0fV , line for any inlet diameter 

would occur in an engine for which the length average internal diameter 
approaches the throat diameter dg. Since dg is related to d^ by 

continuity, this maximum cooling load for the ramjet discussed for the 
supersonic cycle (A^/Ag = 740) is 
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Wing ' * °' 682 a ?’ 8 ^ CIS) 

The total heat added to the system per second is 

E added = u i A l e ( 17 ) 

At a 328,080-foot orbit altitude, with e from. figure 1, this is 

E added * C 18 ) 

Therefore, the desired crude formula for the maximum fraction of heat 
lost by cooling is obtained by combining equations (16) and (l7): 

n v 

Maximum fraction heat lost = — (19) 

R 1 • 

where R^ and are in feet. 

Equation (19) is not a useful rule of thumb for estimating the de- 
parture of the engine from adiabatic, because the ass ump tion of equation 
(16) is probably too severe on the internal design. For example, if the 
engine is 100 feet long and the inlet radius is 20 feet, about 80 per- 
cent of the recombination energy would be unavailable for propulsion be- 
cause of internal cooling. 

An opt im istic value of the cooling load can also be obtained from 
equation (15), which sets a lower limit on the cooling required. By 
assuming that the length average internal diameter is equal to the inlet 
diameter (rather than to the throat diameter), the following equations 
are obtained: 


Pooling = * O'O 49 t-l 8 # 


( 20 ) 


Minimum fraction heat lost = 


0 . 021 - 2 ? 


R: 


1.2 

1 


( 21 ) 


For an engine 100 feet long with an inlet radius of 20 feet, at least 6 
percent of the energy added would be unavailable for propulsion in the 
ramjet example of the cycle analysis. 

The problem of internal friction drag will not be considered, be- 
cause this heat-transfer analysis clearly points out the desirability of 
small ratios of engine length to inlet radius. Designs minimizing in- 
ternal he at losses w ill tend also to reduce the ana logous friction-drag 
loss in the engine. 
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Discussion 

Atmospheric models . - The Rocket Panel model atmosphere (table I) 
was assumed in the preceding analysis. Some consequences of this choice 
will now he discussed. 

Despite the fact that both the cooling and drag problems decrease 
with increasing altitude, a low orbit altitude of about 60 miles was 
chosen for the example, because the kinetics of recombination are sensi- 
tive to total density and to the fraction of atomic oxygen. Since kinet- 
ics are theoretically the deciding factor in the engine efficiency, the 
ambient conditions for a desirable orbit would be high total, density for 
short reaction times and a large fraction of atomic oxygen so that the 
recombination energy would provide appreciable temperature rise in the 
engine. Therefore, this brief discussion of the effect of the assumed 
atmospheric model on the cycle analysis is limited to the recombination 
problem at altitudes where the volume percentage of atomic oxygen is 
appreciable . 

If table II had been used in the cycle analysis, the assumption of 
an orbit at least as high as 360,000 feet would have been necessary in 
order for an appreciable percentage of atcenic oxygen to be present. The 
total density from table II at this altitude is about one-tenth that 
used in the supersonic-cycle analysis. If it Is assumed that the engine 
has about the same ratio of inlet to throat area, then less than 10 per- 
cent of the oxygen can be expected to recombine in a 100-foot nozzle. 
Higher altitudes give less favorable results. Therefore, if a gradual 
transition region in percentage of atomic oxygen like that in table II 
should be confirmed by future research, a ramjet operating only on re- 
combination energy probably would not be feasible. 

The mass fraction of atomic oxygen at 328,080 feet in table III is 
very s imila r to that in table I. However, the ambient pressure and den- 
sity of table III are about 20 percent smaller than the table I values 
at 328,080 feet. An estimate of the effect of these lower ambient con- 
ditions on recombination showed that about 50-percent recombination may 
still be possible in a 100-foot nozzle of an engine with a ratio of in- 
let to throat area of about 740, if kp is 10 13 (g-mole/liter)“ 2 sec -1 

As a final comment, it should be emphasized that the crude kinetic 
analy sis used in this report is justified only by the uncertainty in 
atmospheric composition and experimental rate data. When future research 
clarifies these uncertainties, then the cycle should be analyzed more 
exactly. 

Alternative cycle. - A possible remedy far both the kinetic and 
cooling problems would be to inject a "third-body” gas into the flow at 
the throat. The injected gas would increase the effective density and 
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thus promote recombination in the nozzle. This scheme would detract 

1 from the fuelless feature hut might make the engine operative at highe r 
orbits. Furthermore, a gas-injection cycle would help solve the internal 
cooling problem, because the gas could be used in a regenerative cooling 
system. After cooling the internal surfaces, the hot gas co uld be in- 
jected into the flow to conserve energy lost by the closed-loop radiator 
scheme proposed for the fuelless cycle. 

5 The problem with a gas -injection cycle is how to inject the gas 

2 into the supersonic flow efficiently. If there were shock losses due to 

* injection, the engine efficiency would be lower t itan predicted for the 

supersonic cycle, but the added mass -flow rate in the nozzle would help 
counterbalance this loss . However, the engine would probably act as a 
good brake if a normal shock occurred at the gas -injection station. 


M 
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CONCLUDING REMARKS 

At every turn in this preliminary study, it has been necessary to 
make optimistic assumptions . Seme of these assumptions concern future 
technology, and others concern facts of nature that are not yet known 
with certainty. Furthermore, this report has considered only the steady 
operation of the recombination ramjet in a circular orbit, and the prob- 
lem of starting the supersonic inlet has not been discussed. The thermo- 
dynamic cycle analysis for an ideal engine with supersonic isentropic 
flow throughout and chemical-equilibrium nozzle flow indicated that siz- 
able net thrusts could be realized (fig. 9). It would be possible to 
make allowances for inefficiencies due to internal heat transfer and 
drag, but the major assumption is equilibrium exhaust flow. Taking the 
preliminary results of an AVCO shock-tube experiment as the reccmbinat ion- 
rate coefficient in the nozzle, there appears to be no chance for attain- 
ing chemical equilibrium in the expansion. However, sufficient recombi- 
nation despite departure from, equilibrium may be possible if engine 
lengths the order of 100 feet are used and if appreciable photochemical 
dissociation of air occurs at an altitude as low as 328,080 feet. How- 
ever, drag and heat losses in internal cooling require that the ratio of 
length to inlet radius be kept as small as 5, or preferably smaller, for 
an engine 100 feet long. 

An analysis of a recombination ramjet with normal-shock inlet and 
equilibrium nozzle flow showed that this cycle is not feasible. 

Even though the reccmbinat ion ramjet appears to be a marginal de- 
vice, the problems of this engine do not appear impossible to solve if 
favorable atmospheric composition and reaction-rate coefficients are 
measured in future research. In that case, the fuelless ramjet could 
probably be developed for service where most of the lift is provided by 
the orbiting velocity. 

Lewis Flight Propulsion Laboratory- 

National Advisory Committee for Aeronautics 
Cleveland, Ohio, February 12, 1958 
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APPENDIX A 
SIMBOIB 

A area, sq ft 

atom of "air" 

form-drag coefficient for nonlift elements 
Cp friction-drag coefficient 

Cp heat capacity of air at constant pressure 

concentration, g-mole/liter 
D total drag, lb(F) 

D' drag not due to lift, D - D^, lb(F) 

D l drag due to lift, lb(F) 

d diameter, ft 

E available energy, Btu/sec 

e available chemical-energy density, Btu/cu ft 

F thrust , lb (F ) 

p 

g c conversion factor, lb(M)ft/lb(F)sec a 

H air enthalpy, Btu/lb(M) 

molar air enthalpy, Btu/lb-mole or cal/g-mole 
A heat of reaction, appendix F 
h volumetric air enthalpy, Btu/cu ft 

h„ convective heat-transfer coefficient, Btu/(sec)(sq ft)(°F), 
lb(M)/(sq ft) (sec) in eq. (Dl) 

J mechanical equivalent of heat, ft-lb(F)/Btu 

K chemical-equilibrium constant 

K' ratio of actual concentrations analogous to K, eq. (F9) 

Kn Knuds en number 
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kp recombination -rate coefficient, (g-mole/liter )~^sec”^ 

t 

kp dissociaticn-rate coefficient, (g-mole/liter)~^sec - ^ 

L lift, lb(F) 

££ length of engine, ft 

l mean free path of air 

M Mach number 

Ji molecular weight 

P static pressure, lb(F)/sq ft 

Pr Prandtl number 

q heat -transfer rate, Btu/sec, eq. (ll) 

R radius , ft 

Re Reynolds number 

& universal gas constant 

r recovery factor 

S air entropy, Btu/lb(M)°R 

St Stanton number 

T absolute temperature 

T' lagging temperature (see appendix F) 

t time, sec 

U velocity, ft/sec 

W weight of aircraft, lb(F) 

X any air particle (atom or molecule) 

x length dimension 

y.^ mass fraction of component i, g(i)/g( mixture) 

z ratio of 28.97 to molecular weight of dissociated air 
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T ratio of specific heats 

e emissivity 

t] engine efficiency, eq, (Cl) 

0 cone half -angle of nacelle (fig. 3) 

x thermal conductivity of air 

p air viscosity, lb(M)/(ft)(sec) 

p air density, lb(M)/cu ft, g/cc in appendix P 

a Stefan-Boltzmann constant, 0.1713x10"® Btu/(hr)(sq ft)(9R) 4 

t shear stress, lb(p)/ft 

Subscripts: 

aw adiabatic wall 

C convective 

concentration 

c "combustion" chamber 

e equilibrium 

P frontal projected area of all external surfaces 

f forward reaction 

/ film temperature 

1 any component 

N 2 nitrogen 

n any nozzle station 

0 atomic oxygen 

Og oxygen 

R radiation 

r reverse reaction 
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S solar 

s surface 

T temperature 

w wall 

X any air particle 

y mass fraction 

1 inlet conditions 

2 throat or ambient to flat plate 

3 exhaust -nozzle exit 

oo free stream 

Superscripts : 

~ average 

* reference conditions 
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APPENDIX B 

COMPARISON COP TABULATED PHYSICAL PROPERTIES AND CHEMICAL- 
ENERGY' SOURCES OP UPPER ATMOSPHERE 

From 1946 to 1955, over 200 rockets were fired in upper- atmospheric 
research, and more than half of these parried instruments for measuring 
the pressure and density of the atmosphere. Various authorities in 
upper -atmospheric physics have analyzed the accumulated data and have 
proposed tentative model atmospheres. The following discussion of the 
three atmospheric models presented in the text is included here for con- 
venience. No implication as to the relative merits of the three model 
atmospheres is intended. 


Tabulated Physical Properties 

In 1955, pressure had been observed up to 425,000 feet, and density 
data extended to 725,000 feet. The probable errors in pressure and den- 
sity were less than 10 percent at altitudes below 200,000 feet, but the 
errors increased with increasing altitude up to 200 to 300 percent at 
660,000 feet (ref. 3). 

Below 200,000 feet, the three atmospheric modelB of tables I to III 
agree within 10 percent. Above this altitude, the pressure and density 
of the three tables are difficult to compare, because no general trend 
persists over the entire altitude range. However, in the 300,000-foot 
region with which this report is primarily concerned, table I pressure 
and density are roughly 50 percent greater than table II values and 
about 20 percent greater than table III values. Reference 3 points out 
that a new set of density data was available in 1953 after the Rocket 
Panel had prepared its summary. Measurements of solar soft X-rays be- 
tween 357,000 and 420,000 feet (ref. 14) indicated that the density in 
this region was lower than the average of previous experiments. Table 
II densities over this altitude range are the actual results of. this ex- 
periment, while table III reflects these results with a smaller weighting. 

The temperatures presented in all three tables were calculated from 
the perfect-gas law by assuming a mean molecular weight for air and using 
the given values of pressure and density. The mean molecular weight is 
determined by the chemical composition of air, and at present this com- 
position is still an active research problem. Therefore, the variation 
in tabulated temperatures at very high altitudes is primarily the result 
of the different chemical compositions listed. 
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Tabulated Chemical Composition 

The inert gases , which compose less than 3 percent of the atmosphere, 
are lumped into the tabulated percentage of nitrogen throughout this 
report . 

Atmospheric ozone . - The vertical distribution of ozone from. 30,000 
to 230,000 feet has been measured directly five times with rockets 
equipped with spectrographs (ref. 3). The results are in good agreement 
with the vast amount of data on "total atmospheric ozone" collected by 
ground-level and balloon observations during the past 40 years . The 
averaged results of the rocket observations are given as table XV. 

The averaged concentration data have a probable error of about 20 
percent. Such excellent precision in measuring a "trace" component is 
due to the highly refined ultraviolet absorption technique employed . 

The sharp concentration peak at 100,000 feet is well established, and 
near this altitude ozone contributes s ignif icantly to the available 
chemical energy. 

Gaseous ions ■ - Vertical distribution of electrons has been meas- 
ured with high- alt itude rockets, but physical theories of the ionosphere 
indicate that negative ions are also present (e.g., ref. 15). Therefore, 
the concentration of positive ions cannot be reliably inferred from, elec- 
tron densities. Apparently, ion densities have not been directly ob- 
served in rocket soundings. Furthermore, the energy of ionization might 
not be recoverable in a simple recombination-ramjet cycle. Since the 
ion density is uncertain, and the availability of ionization energy for 
thrust in a heat cycle is also questionable, the contribution of gaseous 
ions to the available chemical energy has been neglected throughout this 
report . 

Atomic oxygen and nitrogen . - The only chemical constituents other 
than molecular oxygen and nitrogen that exist in sizable volume percent- 
ages in the atmosphere are atomic oxygen and nitrogen. 

The Rocket Panel did not seriously propose that the composition of 
the atmosphere corresponds to that shown in table I above 260,000 feet. 

No data of chemical composition were available prior to 1952. However, 
in order to calculate temperature from the perfect-gas law using meas- 
ured pressure and density, the Rocket Panel was forced to assume a mean 
molecular weight for air. Assuming that no dissociation occurred and 
that the molecular weight remained constant, they calculated the left 
temperature column of table I {.JK = 28.97). Believing that these results 
were too high, the panel arbitrarily assumed a uniform dissociation rate 
for oxygen from zero at 262,464 feet (80 km) to complete dissociation at 
393,696 feet (120 km), and similarly far nitrogen from zero dissociation 
at 393,696 feet to completion at 721,776 feet (220 km). The panel 



26 


NACA TN 4267 


presented the second column of temperature (table I, varying ) as the 
"adopted value," emphasizing that they believed the dissociation assump- 
tions were plausible. 

An important set of composition data was obtained in conjunction 
with the total-density measurements from. 357,000 to 420,000 feet in 1953 
(ref. 14). A photon counter sensitive to molecular oxygen was also em- 
ployed in these soundings with Aerobee rockets . The molecular oxygen 
particle density was measured directly throughout this altitude region. 
The percentage composition of atomic and molecular oxygen can be calcu- 
lated from these measured molecular oxygen densities if both the total 
air density and the weight fraction of cxygen to nitrogen are known. 

Table II from 357,000 to 420,000 feet was constructed to correspond both 
to the total air and to the molecular oxygen densities observed in these 
1953 flights, assuming the atmosphere would have been 21 volume-percent 
molecular oxygen if none were dissociated. However, Hulburt (ref. 3) 
points out that the total air density by X-ray measurements of 1953 "gave 
upper-air densities about l/2 to l/3 the values of table I in this alti- 
tude range. One c ann ot decide from the data at hand which is correct." 
For the composition below 357,000 feet, Hulburt used the atomic oxygen 
density results of a photochemical theory that accurately predicted the 
observed ozone density from 60,000 to 250,000 feet. Table IV summarizes 
these theoretical predictions, because atomic oxygen compositions smaller 
than 1 percent are not listed in table II. Above 420,000 feet, the com- 
position given by table H is assumed. The particle densities of these 
constituents above 328,080 feet were calculated for equilibrium with 
gravity. The evidence for the dissociation of nitrogen was mainly 
theoretical. 

As previously mentioned, the total densities in table III are 
slightly greater than the Hulburt -URL values from 357,000 to 420,000 
feet. Therefore, it is not surprising that the ARDC Composition in this 
region differs slightly from that in table II, because higher density 
leads to a different interpretation of the Aerobee soundings . Below 
357,000 feet, a more gradual transition in atomic oxygen was assumed in 
table III than in table II. Above 420,000 feet, table III composition 
parallels that of table II rather closely. 

An interesting discussion of the 1953 Aerobee experiments has been 
published by the Ionosphere Research Group at Pennsylvania State College 
(e . g. , ref. 16). The group presents an analysis that serves as a guide 
to extrapolating the experimental molecular oxygen density to higher 
altitudes. 
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Calculation of Available Chemical Energy 
The following heats of formation were taken from reference 17: 


Particles 

Heat of formation, 
kcal/ g-mole 

at 0° K 

at 298° K 

Ozone, Oj 


34.0 

Atomic oxygen, 0 


59.16 

Atomic nitrogen, N 

112.5 



The heat of formation is the change in enthalpy that would occur if mo- 
lecular oxygen and nitrogen were converted to the listed particles at 
the given temperature . In these three cases , the energy associated with 
the reaction is very large compared with sensible -energy changes. For 
example, the heat of formation of atomic oxygen at 4000° K is 61.161 
kilocalories per gram-mole, or only slightly greater than that listed 
for 298° K. 


Until much better composition data are available, it is permissible 
to neglect the effect of temperature change with altitude and engine com- 
pression and to calculate the "available chemical energy" in the follow- 
ing ma nn er. From table I, at an altitude of 459,312 feet. 


Density, p, lb(M)/cu ft__. 4.747XL0" 0 

Mean molecular weight, , lb/lb-mole "air" 21.21 

Volume -percent atomic oxygen 30.7 

Volume -percent atomic nitrogen 23.1 


The total available chemical energy is equal to 


Ah, 


0-2 


Ah- 




AhQ^ ~ 0 from table XV 

Aho = ^ 5 ‘ - 5 x X 59.16 X (453.6 X 3.968) 

= 7 .32 8X10” 7 Btu/cu ft 

Ah. = 2 - ?u k x 4 ‘ 747X10 — - X 112.5 X (453.6 X 3.968) 

^ 100.0 21.21 


= 10.486X10 -7 Btu/cu ft 




28 


NACA TO 4267 


C 

Therefore, total available chemical energy is 1.78XL0" Btu per cubic 
foot. This value is given in figure 1 at 450,000 feet for the Rocket 
Panel. 
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APPENDIX C 


COMPARISON OF ENERGY AVAILABLE WITH ENERGY REQUIRED FOR FLIGHT 


In this appendix seme relations are developed to examine the flight 
spectrum permitted by the distribution of chemical energy in the atmos- 
phere. Because the engine considered is a ramjet, only the lift-drag 
ratios for supersonic craft will be assumed. Only a rough treatment of 
drag is given here because of the cursory nature of this examination; 
these values are refined in the premising parts of the flight spectrum 
when the ionosphere ramjet is considered in detail in appendix D. The 
object of this analysis is to determine over what range of flight speed 
and altitude a ramjet using the chemical energy of the upper atmosphere 
can produce useful thrust. 

The thrust of the ramjet engine can be written as 


F 


nSJ 

u i 


(ci) 


The lift -drag ratio 
steady flight is 


l/D^ of the wing- fuselage combination in 


If the energy E 
becomes 


_L_ _ W 
D l " F - D' 

is obtained from the intercepted atmosphere. 


(C2) 

it 


E = U 1 A 1 e 

Solving equations (Cl), (C2), and (C3) for e 


yields 


e = 


_W ^L + D 
A 1 L 


A-|_ t]J 


(C3) 


(C4) 


The following breakdown of the nonlift drag terms will 

D' Pi 11 ! p A s p l U l_ ?? P l rM l (% A F 

A ± = °D Aj_ 2g c + °fr A ± 2g c fr 2 \^ fr A 1 


be used: 

$) < C5 > 
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where M-j_ is the flight Mach number. As first approximations , let 


£-1 


and 


As 


= 10 » 


c fr 


> 


(C6) 


J 


A brief survey of current ramjet practice shows that the following 
assumptions are realistic: 


TT- - 100 


lb(F) 

sq ft 


■N 


_L_ 


10 


> 


1-^ 

Rg 


(C7) 




where R is 2 .112X10^ feet, the approximate radius of the earth, and 

p 

the gravitational constant g is 32.2 feet per second . 


Equation (C4) can be simplified by using equations (C5) to (C7) as 
follows : 


e = 


T]J 


10 1 - 


2 

Rg 


|+ 5P-J-; 


fr 


(C8) 


Equation (C8) gives the minimum energy density_necessary to sustain 
flight in terms of the flight conditions and C^ r . 


The flow conditions in which the ramjet wing-fuselage will be fly- 
ing must first be determined before the friction-drag coefficient is 
evaluated. Although the boundaries of slip and continuum flows are not 
sharply defined, the definitions of reference 18 will be adopted here. 
That is, the limit of continuum flow occurs where the Knuds en number 
(Kn = if <£) is less than 0.001, and slip flow exists for Kh from 
0.001 to 2. If the characteristic length of the ramjet is 10 feet and 
the total length at least 100 feet, then the flow will be slip or con- 
tinuum below about 425,000 feet. Furthermore, in the slip-flow region, 
Cfr £° r continuum flow will probably be within 10 to 20 percent of the 

correct value (refs. 19 and 20) . Therefore, C^.. will be calculated 
from an empirical method that is applicable to continuum laminar flow. 
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The length average friction-drag coefficient for laminar flow over 
a flat plate is 



(C9) 


Eckert (refs. 6 and 7} has shown that, by evaluating the thermodynamic 
and physical properties of air at a defined reference temperature, equa- 
tion (C9) can be -used for a wide range of wall temperature and Mach num- 
ber conditions . This calculation method is discussed in detail in appen- 
dix Dj for approximate results, the reference -temperature method is con- 
venient. The reference temperature for evaluating p and p is defined 
as 

T* = 0.5(T w + Tj) + 0.22(T aw - Tj (CIO) 

where 

V - + Hr rM ?) 


and 


T. 


r = 


aw 


- T, 


Ui 


2Jg c c p 


It is convenient in this approximate analysis to use free-stream values 
of p and p, designated by the subscript 1. Therefore, the effect of 
temperature on air properties is approximated by 


and 



so that 



(Cll) 


Before the friction-drag coefficient determined by Eckert’s method 
can be applied to this analysis, the definition for C fr by Eckert must 



32 


M.CA TH 4267 


be reconciled with the definition used in equation (C5). That is, the 
reference-temperature method defines the average shearing stress i as 


or 


* ■ I 


-In tt 2 

T “ 2 C fr p i u l 



(C12) 


However, in equation (C5) the implicit definition is 


\ 


(C13) 


The relation between the Eckert C fr in equation (C12) and the C fr 


i simply 


- f T l\ 


M?) 

= C fr 

eq. (C12) 


eq. (C5) 


(C14) 


Combining equations (C9) and (Cll) yields 

eq. (C12) 


'fr 


= 1.528 f T l^ 

V 


(C15) 


Finally, the desired coefficient for equation (C5) results from, substi- 
tuting equation (C15) into (C14): 


'tx 


eq. (C5) 


, . 0.2 

_ 1.528 f ^l\ 

’ . OOE\t*/ 


(C16) 


IT 

p l 


Consider the friction-drag coefficient for two flight-speed condi- 
tions. Low supersonic Mach number: 


T- 

T' 


1 ■ 


(C17) 


High supersonic Mach number: For orbital Mach number of about 23, 

T aw/ T l = 1 + 0-2(550); and, assuming a material limit for the wall 
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in 

t 


B 


temperature, Ty/T-^ - 5. Thus, equation (CIO) can "be evaluated as 


f~ = 0.5(5 + 1) + 0.22(110) = 27.2 


and 


✓ m x 0 - 2 

(?) 


1 

2 


(C 18 ) 


Furthermore, it is assumed that 

PlU x x 

*1 


= ICf 


(C19) 


Then, the C fr from equations (C17), (C19), and (C16) and from 
(C18), (C19), and (C16) are 


C fr = 1.328X10 


.-2 


for low Mach number 


= 0.664XL0 2 for high Mach number 
Finally, substitution of equation (C20) into equation (C8) gives 


(C20) 


= ^ 10 ^1 - ^1.328XL0 -2 ccr 0.664X10~ 2 ^ 


(C21) 


At low supersonic velocities and at high altitudes, the first term will 
dominate . Taking q = 100 percent, 


e = -^| = 1.29XL0 -2 


Comparison of this value with figure 1 shows that the atmosphere cannot 
provide the required energy anywhere. 

For orbiting conditions, the first term of equation (C2l) vanishes, 
and taking y = 1.4 and ~ 23, the equation becomes 

e = (3.27XL0" 2 )P 1 
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Using table I densities and pressures together with figure 1 , the follow- 
ing values are obtained: 


Altitude, 

ft 

Required 

e 

Available e 

Low 

High 

100,000 

7 .3X10" 1 

4.0XL0 -5 

4.0XL0" 5 

300,000 

1.4XL0" 4 

4. 0X10" 6 

9. 0X10" 8 

350,000 

1.7XL0' 5 

3. 0X10" 6 

2. 0X10" 5 

400,000 

2. 9X10" 6 

1. 0X10" 6 

5.0XL0" 6 

500,000 

2. 9X10" 7 

1.0XL0' 7 

9. 0X10" 7 

600,000 

4. 9X10" 8 

5. 0X10" 8 

2. 0X10" 7 


Thus, above 300,000 feet there is the possibility of obtaining thrust to 
overcome drag with an orbital ramjet that obtains its energy from, the 
energized air it flies through. These values are for a Reynolds number 
4 

of 10 . The influence of Reynolds number is discussed in the text under 
THE IONOSPHERE RAMJET. 
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APPENDIX D 

HEAT TRANSFER AND FRICTICN DRAG 
By James F. Schmidt 


Eckert's Laminar Heat-Transfer Theory for a Realistic Gas 

Eckert's laminar flat-plate solution for the heat-transfer rate and 
skin-friction coefficient is applied herein to the three configurations 
shown in figure 3, a truncated cone, a circular cylinder, and a reversed 
truncated cone. Since large temperature variations are evident through- 
out the boundary layer at high speeds, Eckert's reference-enthalpy method 
is used rather than his reference- temperature method (refs. 6 and 7). 

This reference -enthalpy method accounts for the varying specific heat in 
boundary layers with large temperature variations. 


The predominant reason for using Eckert's theory is that his approach 
is a direct calculation method that gives results comparable to those of 
methods that assume equilibrium dissociation of air (ref. 2l) . Reference 
21 shows that the skin friction and heat transfer are relatively unaf- 
fected by dissociation if the surface temperature does not exceed 2500° R. 
In recent years some thought has been given to the phenomena of diffusion 
of dissociated particles. Some investigators have attempted to estimate 
the diffusion rates that could be used in heat-transfer calculations. 
However, sufficient information on diffusion rates is not available for 
such calculations at the present time. 


All the heat -transfer calculations herein are based on the atmos- 
pheric properties of table I. The heat -transfer rate is defined as 


1 = VEaw “ 


(Dl) 


where is the enthalpy due to wall temperature, and is the re- 

covery enthalpy, defined as 


Saw = % + r 



(D2) 


where Hg is the enthalpy due to the temperature of the gas surrounding 
the flat plate, and Ug is the velocity of the gas relative to the flat 
plate. The recovery factor r Is obtained frcm the following equation: 
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■where 



(D4) 


A value of Pr* = 0.71 Is assumed in calculations for all configurations 
(ref. 7, fig. 3). 


The reference enthalpy is calculated from, the following semiempir- 
ical equation: 

H* = 0.5(11^ +1^)+ 0.22(^ w - Hg) (D5) 

The reference temperature is determined from the reference enthalpy. 
In the high-temperature range, a linear extrapolation of temperature 
against enthalpy from the gas tables of reference 22 is used. This lin- 
ear extrapolation was checked against the tables of reference 23 for un- 
dissociated air up to 9000° R. Even at 9000° R the accuracy of the lin- 
ear extrapolation is within 5 percent. 


The skin-friction coefficient is given by the following equation: 


o 0.664 



(D6) 


where 

[i* viscosity evaluated at reference temperature T*j formula for vis- 
cosity in UBS tables (ref. 23, p. 10) is -used for all 
calculations 

p* density evaluated at reference temperature T* 
x length dimension of flat plate 

The Stanton number is related to the friction coefficient by 


St = (Pr*)" 2 / 3 (D7) 

u 

and 

p*u 2 


St = 


(D8) 
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Equilibrium External Wall Temperature 

The equilibrium wall temperature is determined from, t he following 
energy balance equation: 

SC + % = C D9 ) 


where 

< 3 g convective heat-transfer rate to surface 
qg solar heat-transfer rate to surface 
q R radiation heat-transfer rate from surface 

Convective heat -transfer rate . - Eckert's flat-plate solution is 
vised for calculating the convective heat-transfer rates. The wall tem- 
perature has only a slight effect on the convective heat-transfer rate, 
pri mar ily because the recovery temperature is extremely large in compar- 
ison with the relatively low wall temperature. 

Solar heat-transfer rate . - The ramjet vehicle is assumed to be 
orbiting but not rotating, so that the sun radiates heat to only half 
the vehicle surface area at all times . In addition, the max imum heat 
energy is radiated to the ramjet by assuming that the sun's rays are 
always normal to the exposed surface area of the ramjet. Since the 
highest equilibrium wall temperature tends to produce the maximum ex- 
ternal drag, the maximum heat energy transferred to the vehicle repre- 
sents the worst flight condition. 

The solar heat -transfer rates far the sunny and dark side of the 
vehicle are assumed to be 0.118 and 0 Btu/(sq ft) (sec), respectively. 

All the solar energy radiated to the vehicle surface is assumed to be 
completely absorbed. This assumption of complete absorbtion of the max- 
imum amount of radiated solar heat energy predicts the largest external 
drag for all configurations . 

Radiant heat-transfer rate . - The heat energy transferred away from 
the vehicle is assumed to be radiating to dark space. Therefore, this 
radiant heat-transfer rate is defined by the following equation: 

OR = aelj (DIO) 

where the Stefan-Boltzmann constant a is 0.17 13X10" 8 Btu/(hr)(sq ft) 
(PR)*, and the emissivity e is taken as 1.0. 
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Equilibrium -wall temperature against altitude . - A balance of energy 
between the convective plus the solar heat flux and the net radiation 
heat flux determines an equilibrium wall temperature (eq. (D9)). An 
equilibrium wall temperature was determined for the case in which n.n 
the energy of the atomic oxygen of the ambient air is assumed to be added 
to the free -stream energy. It is assumed that this condition would pro- 
duce the highest equilibrium wall temperature for any given configura- 
tion . This high-energy condition was applied to configuration I (fig. 
3(a)) with the following dimensions: a truncated ccne with a 2° angle 

of inclination and a total length of 100 feet. A max imum equilibrium 
wall temperature of slightly less than 1000° R was calculated far con- 
figuration I. The combination of these dimensions for configuration I 
and the high-energy condition is believed to give the highest possible 
equilibrium wall temperature. 

In order to ensure more than marginal safety, a wall temperature of 
1000° R is assumed for all drag calculations. Since configuration III 
appears to give the lowest external drag, an equilibrium wall temperature 
was calculated for this configuration at 328,000 feet. This calculated 
temperature was found to he more than 300° R lower than the assumed value 
of 1000° R. Therefore, it is evident that the assumed w all temperature 
of 1000° R is conservative. 

Also, the effect of altitude on wall temperature is negligible for 
altitudes up to 700,000 feet (fig. ll). 


External Surface Drag Per Inlet Area 

Surface drag per inlet area for configuration I in continuum flow . - 
An oblique shock wave is produced at the inlet of the truncated cone 
(fig. 3) . The flow properties behind the shock are found from oblique - 
shock-wave theory (ref. 24) for a perfect gas using free-stream condi- 
tions. The skin-friction coefficient is calculated as a function of the 
flat-plate length x by use of Eckert's flat -plate solution with the 
flow properties behind the oblique shock wave. The truncated cone is 
treated as a flat plate with infinitesimal width at an angle of attack. 
Expressing the surface area of the cone as a function of x, the surface 
drag is obtained by integrating the product of the local friction drag 
and the unit surface area over the surface length of the cone. This 
drag is divided by the inlet area in order to give the external surface 
drag per inlet area: 



(Dll) 


OfiAt 
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where 


dA- s = f (x) = 2 jc(R}_ + x tan 0)dx 

c fr P*u| 

V " 2 S C 

and 

C fr = f(x) 


Therefore, the surface drag per inlet area is 


D 

A i p i 



°fr< E l 


+ _x tan 0)dx 


(D12) 

(D13) 

(D14) 


(D15) 


The pressure drag has negligible effect on the total surface drag for 
this configuration with small cone angles . All drag calculations are 
based on the atmospheric properties of table I. 


Surface drag per inlet area for configuration II in continuum, flow . ■ 
Since configuration II, a circular cylinder, as shown in figure 3 is 
orientated at'0° to the free-stream flow direction, no shock waves will 
occur on the cylinder. Hence, free-stream flow properties are used in 
Eckert's flat-plate solution for the skin-friction coefficient. 

The circular cylinder is treated in the same manner as the truncated 
cone (configuration I) for obtaining the surface drag per inlet area 
(eg,. (Dll)), where 


dAg = f(x) = 2 ^ 2 . dx 
and equations (D13) and (D14) apply. 

Hence, the surface drag per inlet area becomes 

P pVe r 

¥i " g c R iPi J. 


dx 


(D16) 


(D17) 


Surface drag per inlet area for configuration III in continuum 
flow . - A Prandtl-Meyer expansion is produced at the inlet of the re- 
versed truncated cone (fig. 3). Since the expansion angle is small 
(0 = -4.3°, -8.617°), the two-dimensional Prandtl-Meyer expansion is 
assumed to be a good approximation for the flow around the corner of 
the cone. 
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Using the properties behind the Prandtl-Meyer expansion in Eckert's 
flat-plate solution, the skin-friction coefficient is calculated as a 
function of the flat-plate length x. The reversed truncated cone is 
also treated in the same manner as the truncated cone for calculating 
the surface drag per inlet area (eqs. (Dll) to (D15)). The pressure drag 
for this configuration is neglected, because it actually produces a neg- 
ligible force in the thrust direction. 


Maximum Internal Heat-Transfer Rate 


The maximum internal heat-transfer rate is encountered at the throat 
of the ramjet nozzle. The internal heat -transfer calculation was made 
for configuration III with the following geometric dimensions: 


9 = -8.616° 
& = 100 ft 



d ] _ = 40 ft 
A 1 

£ - = 740.65 
dg = 0.465 ft 

The heat-transfer rate is defined by the following equation: 


1 - 


(D13) 


The following values were assumed: Allowable wall temperature with max- 

imum internal cooling T w , 1660° R, and adiabatic temperature of the 

dissociated air T &w , 9850° R. 


The heat-transfer coefficient Ilq (ref. 13) is determined from 

.0.8,„ \ 0 . 33 


h0 " \ <°- 023) 



(D19) 


/ 


where the film temperature T y is 3311° R , and 

X y = 174. 5X10” 7 Btu/(ft)(sec)(°R) 
p 2 = 410X10" 7 lb/cu ft 





4=780 


NACA 131 4267 


41 


U 2 = 24,743 ft /sec 
iy= 419XL0 -7 Ib/(sec)(ft) 
(GpMA ) y = Pr = 0.71 
The resulting heat-transfer rate Is 

JL = 10.82 Btu/(sec)(sq. ft) 


& 
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APPENDIX E 


NORMAL-SHOCK ORBITING RAMJET CYCLE 

The all-supersonic ramjet cycle discussed in the text gives ideal 
engine efficiencies up to 85 percent. However, since the problem of 
starting such a ramjet might be great, it is natural to investigate cy- 
cles involving inlet shocks. The extreme case of a normal- shock inlet 
will be examined here. The lower engine velocities would give longer, 
more favorable reaction times j and, as outlined in appendix F, the re- 
combination kinetics would be a critical problem. 

For this example, the flight altitude and velocity are identical 
with those of the all- supersonic cycle. The problem is twofold: (l) 

The flow conditions after shock must be determined, and (2) the cycle 
must be completed by expanding the air to ambient pressure and obtain- 
able thrust must be calculated. 

The flow immediately following the normal shock will not be in chem- 
ical equilibrium, because even dissociation kinetics are relatively slow 
at these rarif ied gas conditions . However, the air after shock will ap- 
proach equilibrium if isentropic constant-area flow sufficient for reac- 
tion is assumed to exist downstream. This assumption is made in the fol- 
lowing analysis. 

The following four equations .must be satisfied at the inlet: 
Continuity: 


PooUoo - P]_U]_ 

Conservation of momentum: 

p i - “ 7T (Pj£ " PiU?) 

°c 


Conservation of energy: 

2g c J( Hl - Hj = U* - U* 


Equation of state: 


P 1 = ^ z lPi T l 


(El) 


(E2) 


(E3) 


(E4) 


The subscript « refers to the free-stream conditions, and 1 is the 
inlet or after-shock flow. The solution of these equations is straight- 
forward, but it involves a trial-and-error procedure and the use of 
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o 

oo 

F- 


charts of thermodynamic properties of dissociated air. The known free- 
stream. conditions are identical with those for the all- supersonic ramjet 
cycle used in the text: 

U*, = 26,050 ft/sec 
p ro = 5.375XL0" 8 lb(M)/cu ft 
P*, = 1.24XL0 -3 lb(F)/sq ft 
= 984.8 Btu/lb(M) 


z = 1.103 

CX> 

Briefly, the calculation procedure is to use an assumed in a 

combination of equations (El) and (E2) and also in equation (E3). Fur- 
thermore, it is assumed that the air after shock is in chemical equilib- 
rium, so that the calculated and H-j_ dete rmin e T]_ and Zj_ on a 

Mollier chart for air (ref. 8). Equation (E4) with the determined val- 
ues of z^, T^, and P^ yields the corresponding p-, . Finally, this 
procedure is repeated until the continuity equation (El) is satisfied 
by the assumed and the calculated p-j_. In this manner, the flow 

conditions after shock were found to be 


U ± = 1195 ft/sec 
Pj_ = 1.13 lb(F)/sq ft 
p x = 1.17X10“ 6 lb (M) / cu ft 
T-l = 9630° R 

Zj_ = 1.888 

E 1 = 14,400 Btu/lb(M) 


*^ 1 S 1 

st 


= 65.55 


The r ema inder of the cycle is easily determined on the Mollier chart by 
following the isentropic curve 31 to ambient pressure or P ro . 

VhmigrVi the chart used did not extend over the entire pressure range, a 
crude extrapolation is sufficient; this will be evident after the thrust 
calculation . 
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The nozzle exhaust properties by extrapolation are 

T 3 = 3650(1.8) °R 
H 3 = 7250 Btu/lb(M) 
z 3 = 1.45 

From the energy equation and continuity, 

U 3 = 18,950 ft/sec 
A* 

r - 

00 


The thrust from this normal -shock- inlet cycle was calculated from 


F 


Poo^a 

g„ 


J+*s(r)- 

\ 00 / 


O) 


» -10 lb(F)/sq ft 

This large drag clearly reflects not only the inefficiency of a shock 
cycle but also the fact that the great dissociation at the extreme inlet 
temperature cannot be recovered in the exhaust nozzle to any appreciable 
extent . - 
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APPENDIX F 


CHEMICAL KINETICS QE EEC CMBLJATICSN 

The first two sections of this appendix apply the general method of 
reference 10 to the recombination of atcmic oxygen in an adiabatic noz- 
zle flow. An equation that is used as a criterion for near-equilibrium 
flow is obtained first, and then a similar relation is obtained for near- 
frozen expansion. In the final section, these limiting cases are con- 
sidered in a numerical example to point out the necessary kinetic condi- 
tions on the thermodynamic cycle assumptions. 


Adiabatic Expansion with Near-Equilibrium Flow 

The intuitive meaning attached to "near-equilibrium flow" is that 
the difference between the actual weight fractions of a given chemical 
component y.^ and the equilibrium weight fractions for the local condi- 
tions of temperature and pressure y^ e is not large. Before stating 

the precise definition, it will be necessary to review some background 
thermodynamic and kinetic equations. Then, a -useful criterion for near- 
equilibrium. flow will be reduced from a mathematical formulation of the 
preceding qualitative statement . 

For this analysis, let the actual molecular reaction be 

0 + 0 + X ^ 0 2 +X (Fl) 

k r 

The rate of the forward recomb in ation reaction is given by 

where <£ ^ is the concentration of species i (g-mole(i)/cc) and the 

subscript 0 refers to atomic oxygen and X to all "third bodies," 
taken here as the total molar density. This expression assumes iso- 
thermal constant- volume reaction. Similarly, the dissociation rate is 

(K5) 

where the subscript 0g refers to molecular oxygen. 
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The net production rate of atomic and molecular oxygen in the re- 
versible reaction equation (Fl) can be written in terms of equations 
(F2) and (F3) as 

1 d ^ 0 dSf 0 

dt dt“ - ‘ " CF4) 

The factor l/z is determined by the stoichiometry of the reaction; two 
moles of atomic oxygen react for every mole of molecular oxygen formed. 


Later in this analysis, it will be convenient to work with the con- 
tinuity equation in terms of weight fractions rather than molar concen- 
trations . Therefore, noting that 


y ± p 

jT^ 


(F5) 


where 

y^ weight fraction species i, g(i)/g (mixture) 
p total mass density, g(mixture)/cc 

molecular weight of i, g(i)/g-mole(i) 
equation (F4) is rewritten: 


d'tfg 

u 2 

dt 


1^0 
2 dt 



At this point, the thermodynamic chemical equilibrium constant at 
the local temperature and total pressure can be introduced. Defining 
the equilibrium constant in terms of equilibrium concentrations in the 
usual manner and then using equation (F5) to obtain its relation to an 
equilibrium constant in weight fractions give 
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The subscript e denotes that equilibrium values must be used. Further- 
more, at equilibrium, the net rate of production of all chemical species 
is zero. Thus, 

' e (F8) 


K - k f - 

**■ " 35 ■ 


(e’oj 


Equation (F6) can be written in terms of Ef and by using equa- 

tion (F8) . However, before doing this , it is convenient to define a 
term analogous in form to the equilibrium constant but expressed in terms 
of the actual weight fractions at the local temperature and pressure: 

K * = — -£■ 


V 






Q T 

K <r = ^Tp^ 


> 


(F9) 


Substituting equations (F8) and (F9), 




d^ 0 . 


dt 


l d ^0 
2 dt 


= kf 




2 

y^x 


(■•«) 


(fio) 


Note that the subscript has been dropped from K'/K, because from equa- 


tions (F7) and (F9) it is apparent that 


K. 




K y 


K 


% % 


K' 

K* 


Reference 10 shows that a one -dimensional-flow continuity equation 
for negligible diffusion currents is 


d m y02 

Dt 


_Dlny 0 1 d ^0 2 ^ x 

2 Dt <<f 0 dt " <# Q dt 

2 


(Fll) 


where the Euler total derivative for steady x-directicn flow is 
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A coordinate transformation on .the left side of equation (Fll) introduces 
the time rate of change of temperature T: 


M In y 0g ^ d In y Q \ m 1 d ^0 2 ^ ± d«r 0 

\ dT " dT J Dt = <T 0 dt ~ <# a dt 

u ^ 


(F12) 


The desired equation for developing criteria, for near-equilibrium, flow 
and for near-frozen flow is obtained by combining equations (F9), (FIO), 
and (F12): 



The statement at the beginning of this analysis indicated that in 
near-equilibrium flofr the difference between y^ and y^ e is not 

large. Following Penner, the quantity Ky can be considered to be an 

equilibrium constant at a temperature T' that is slightly different 
from the actual local temperature T. The precise definition for near- 
equilibrium flow states that the temperature difference "T 1 - T is so 
small that it is sufficient to retain only the first two terms in the 
Taylor series expansion of Ky = Ky(T’) about the temperature T " (ref. 

10). That is, for near-equilibrium flow, 


Ity = %(*’) * %(T) + Ky(T) 
Equation (F14) is easily rearranged to 


~d In K„(T)“ 
dT 


(T' - T) 


(F14) 


Ky fd In Ky(T)“ 

% = ’ L ® : 


T 


(T l - T) 


(F15) 


{< d In kA 

Equation (F15 ) , along with the approximation that ( — J 

d In K V (T)1 

} is now substituted into equation (F13), which simplifies 

__ dT J T 


to the following equation in terms of molar concentrations: 


pi _ p at 


(•S) 


k f^x l 49 ’o + 


^ o 


l-l 


(F16) 
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Equation (F16) relates the temperature lag associated with the composi- 
tion of the near-equilibrium reacting gas mixture to the temperature rate 
of change with time at a given point in the exhaust and a reaction t ime 
(in parentheses). 


Adiabatic Expansion with Near-Frozen Flow 


Penner (ref. 10) defines near-frozen flow as flow for which the 
extent of chemical reaction is so small that the difference T c - T' 

always remains small. Here, T c is the temperature of the gas at the 
start of expans ion j and, as in near- equilibrium flow, T' is a temper- 
ature for which Ky(T') = Ky. Therefore, in this case the criterion is 

t d In K V (T)1 

C* 0 - I') (KL7) 


Equation (F13) can be written in terms of molar concentrations as equa- 


tion (F18) by solving equation (F17) for 
approximation that 




and using the 



ThiR equation can be reduced by noting that, in the numerical example 

Ky(T ) 

being considered, <<1 and that T c - T' is also small, so that 

equation (F18) simplifies to the following: 
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Chemical thermodynamics can he used to express the temperature - 
dependence of the weight -fraction equilibrium constant Ky in terras of 
the heat of the oxygen recombination reaction The exact form of 

this relation is complicated, but the following equation is a good 
approximation: 


d 111 % 

W ^T 2 

Therefore, equation (F18) can be written as 


(F20) 


T„ 


T' 


T c - T 



kf ^ 




(F21) 


Equation (F2l) is the desired expression for characterizing near-frozen- 
composition expansion. 


* 




s 


Numerical Example 

Before applying equations (P16) and (F2l) to the nozzle example of 
the cycle analysis, two basic assumptions should be pointed out. The 
most important is the assumption that the reacticn-rate coefficient kf 
is not temperature -dependent . In general, chemical reaction rates are 
very temperature-sensitive. However, for the recombination reaction 
given by equation (Fl) (forward reaction), only an algebraic temperature- 
dependence should be expected if oxygen recombination follows the trend 
observed in iodine atomic recombination and hydrogen atomic recombina- 
tion. In the latter two reactions, near-zero activation energies have 
been observed. For lack of sufficient data the k^ has been assumed to 
be temperature-independent . Secondly, the total derivative DT/Dt is 
assumed to be constant throughout the nozzle; this is an adequate approx- 
imation for the limiting cases being considered. 

Near-equilibrium expansion . - The temperature, pressure, and compo- 
sition conditions in the equilibrium expansion cycle calculation are 
used here. Specifically, 



assuming that 


Se = 100 feet 
U = 25,000 ft /sec 
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g 

c~- 

-=H 


- = 5.5XL0 5 °K/sec 

The molar concentrations of 0, 0g, and X (total) are calculable from the 
equilibrium nozzle conditions given by figures 5 and 6. By assigning a 
maximum, temperature lag (T 1 - T) , equation (F16) was solved for minimum 
values of kf . The results of this .calculation using T' - T = 20° K 
are given in table V. This temperature lag is somewhat arbitrary, but 
an examination showed that 20° K is small enough to assure the rapid con- 
vergence of the Taylor series (eq. (F14)) for the oxygen recombination 
reaction. 


Near-frozen expansion . - The calculation procedure for near-frozen 
expansion is identical to that for near -equilibrium expansion, except 

( T - T'\ 

■fpS — J i 10“ 2 was used in equation (F2l) to calculate the maximum 

kf possible before the flow deviates from frozen expansion. The heat of 
reaction given in appendix B (AJf’ = 2(59,160) cal/mole 0g) and T c = Tg = 

2763° K are the additional information required. The results of these 
calculations are also given in the text. The reduced temperature ratio 

of 10” 2 is sufficient for the convergence of the Taylor series (eq. 

(F17 ) ) to within 10 percent in two terms . 

An apparent inconsistency in the derivation of equation (F2l) and 
its application in this example should be pointed out and explained. The 
derivation ass ume s that the gas at the start of expansion is in chemical 
equilibrium; this justifies the use of equation (F20). The results of 
this crude analysis, as discussed in the text, show that equilibrium 
flow may never be attained in the example nozzle. However, from figure 
6, it is apparent that in this cycle the difference between the ambient 
frozen composition and the equilibrium composition at Tg and Pg is 

very slight, being off only by the presence of 2 percent nitric oxide at 
equilibrium. Therefore, the use of equation (F21) appears justifiable 
in this example. 


Pe rcent reci’ttnhin^a . - The molar concentrations of atomic oxygen cor- 
responding to the necessary a given kj* and flow condition can 

easily be related to the degree of recombination in the nozzle. Since 
the pressure aryl temperature change markedly in the nozzle, the percent- 
age recombination is defined here in mass units: 


Percent recombined = 100 


( y 0^2 " ^ y 0^n 

^ 0*2 


(F22) 
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This equation can be rewritten in terms of the molar concentrations at 
stations 2 and any other nozzle condition n: 


Percent recombined = 100 


2 (^O^n 2 

^n "^0^2 l ^Vn 


(F23) 


The nozzle conditions summarized in figures 5 and 6 were used in equa- 
tion (F23) to calculate figure 10. 
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TABLS I. - ROCKET RESEARCH PAMKLl KDOPKD AOTOSmWIO MIA (RSI . 2) 
[Converted to engineering unit*.] 


Altitude 

(above 

■ •& level). 


Pbjalwa pr°p»rtl*» 



Owileal oa^oiition® 



Te«p, tor 

Preisun, 




VdIujk; 

Volume 

YoIum 

Volume 

Weight 

ft 

00 

3- 28.87, 
°F 

lb (F)/»q ft 



percent 

02 

ptrc»ct 

0 , 

peraont 

*2 

peraont 

q 

ratio. 

Oxygon 

Nitrogen 

5,969 

64.1 

1840 

B.Bsaxia-® 

04 .1 

£8.97 

21 

0 

79 

0 

0.5020 

15,125 

31,0 

1518 

B.01BX1Q-X 

31,0 

28.87 

21 

0 

79 

0 

.5020 

28,246 

-18,7 

777.7 

3.507X10-2 

-18.7 

23.97 

£1 

0 

79 

0 

,5020 

39,370 

-34 .6 

431.4 

2.0*8X10-2 

-84.8 

28.97 

21 

0 

i» 

0 

•3020 

32,493 

-85.5 

227.4 

1.139X10" 2 

-05.3 

28.97 

21 

0 

79 

0 

.5020 

86,818 

-78.7 

118.0 

6. 626X10-® 

-76.7 

£8.97 

21 

0 

79 

0 

,5020 

78,740 

-62.1 

63.66 

3 .002X10"® 

-02.1 

26.97 

21 

0 

79 

0 

,5020 

91,865 

-48,4 

34.82 

1,580X10“® 

-48.4 

20,97 

21 

0 

79 

0 

.3020 

104,980 

-35,5 

IB .46 

8.898X10-* 

-35.5 

28.97 

21 

0 

70 

0 

.3020 

118,110 

-15.1 

11,11 

4 .671X10-* 

-15.1 

£8.97 

5LL 

0 

79 

0 

,3020 

151,232 

12.9 

B .559 

2.B02X10-* 

12.8 

£8.97 

£1 

0 

79 

0 

,3020 

144,558 

27.2 

5.961 

1 .525X10 "* 

27.2 

20,97 

21 

0 

79 

0 

,3020 

157,480 

29.2 

2.408 

9.255x10-8 

29.2 

28,97 

£1 

0 

79 

0 

,5090 

170,805 

25,7 

1.468 

6 .867X10“® 

25,7 

20.97 

21 

0 

79 

0 

.3020 

183,728 

14.9 

.689 

3.511X10- 5 

14.9 

28,97 

£1 

0 

79 

0 

.3020 

196,348 

-4.7 

8.50x10-} 

2.180)0,0-? 

-4.7 

28.97 

21 

0 

79 

0 

.3020 

209,971 

-29,6 

3.07X10"* 

1 .538X10"? 

-29. 6 

£8.97 

£1 

0 

79 

0 

.3020 

223,094 

-se.o 

1,72x10"; 

3.024X10“” 

-68.0 

20.97 

21 

0 

79 

0 

.3020 

258,217 

-75.2 

9,55X10-2 

4.543XX0 - ® 

-76.2 

28,97 

SI 

0 

79 

0 

.5020 

£49,540 

-85.6 

4.95X10" 2 

2.400X10"** 

-85.5 

28.97 

21 

0 

79 

0 

.3020 

262,464 

-90.7 

2.BSK10-* 

1. 518X10 -® 

-90,7 

28.97 

£1 

0 

79 

0 

.3080 

278,868 

-85,3 

1.15X10, 

0.7B1XLO-; 

-93.2 

28.23 

17.9 

5.1 

77 

0 

.3020 

298,272 

-69.1 

5.51X10" 5 

2.540X10"' 

-SB. 6 

27.52 

15.0 

10,0 

75 

0 

.5020 

Ml, 878 

-60.2 

£.51X10-5 

1.151X10-7 

-80.1 

26.66 

12.8 

14.6 

73.2 

0 

.3020 

320,080 

-27.7 

1.24X10"® 

B.375X10" 8 

-88.6 

26.22 

9.5 

19.0 

71.3 

0 

.3020 

344,484 

-7.4 

e.ssxio-* 

«, 611X10"? 

-54.2 

25.61 

7.0 

£3,2 

69.8 

0 

.3020 

360,888 

26.3 

3.30x10"* 

1. 202X10“° 

-59.8 

20.03 

4.5 

27.2 

88.3 

0 

.3020 

377,292 

BO. 3 

1.58X10"* 

8.822X10“® 

-8.8 

24.48 

2.2 

S1.0 

88. B 

0 

.3020 

393,696 

154.5 

1.11X10-* 

5. 518X10-® 

31.3 

25>B6 

0 

54.7 

85.3 

0 

.3020 

410,100 

168. S 

6. 92X10 - b 

1. 997X10-® 

58.5 

25.19 

0 

35.6 

SO.l 

6.5 

.3020 

426,504 

242,3 

“4 ,47KL0'g 

1.190X10'® 

85.5 

22.50 

0 

32,6 

56.2 

12.2 

.3020 

442,906 

294.5 

2 ,90X10"? 

7.388x10’}" 

107.5 

21.85 

0 

51.8 

50 .5 

17.9 

.3080 

439,312 

344.9 

2 .06X10"? 

4.747x10-}? 

12B.4 

21.21 

0 

50.7 

4B.2 

£5.1 

.5020 

476,716 

3BG,S 

1 .45X10"? 

5.158X10-}? 

148.5 

20. &e 

0 

29,8 

42,1 

28.1 

.3020 

492,120 

446.7 

1.05XLO" 3 

2.126X10" 10 

151.4 ■ 

20 .06 

0 

29,0 

56. 2 

52.8 

.5020 

508,624 

498,1 

7,80)00"® 

1.467X10-}? 

183.0 

19.49 

0 

2B.5 

34.5 

37.2 

.5020 

524,926 

546.5 

5. 57X10 J 

1.034X10 

6.412XL0“ j “ l 

202,7 

19.05 

0 

27.5 

81.1 

41.4 

.5020 

857,736 

653.9 

5 . 21X10 

£36,3 

19.10 

0 

£6.2 

£4.0 

49.2 

.3020 

990,344 

758.7 

1. 95X10"® 

5.002X10-}} 

288.4 

17.28 

0 

24. B 

18.8 

58.3 

.3020 

623,332 

665.1 

1 ,24X10-1 

1.751X10-}} 

295. S 

IS. BO 

0 

25.8 

18.3 

62.7 

.5020 

056,180 

966.8 

g.iixio-Z 

1.066X10"}; 

M8.1 

16.79 

0 

22. Q 

8.8 

86.8 

.3020 

883,968 

1069.9 

5 .4BXI0-J 
3. B2X10"' 

6.7&OXLO"} 2 

4.S79XLO - - 1 - 2 

540.2 

16.15 

0 

21,9 

4.1 

74.0 

.3020 

721,776 

1172.2 

300.2 

14.55 

0 

£1.0 

0 

79.0 

.5020 


■Aaeuned from theoretical tmloul.tloos «bov» £62,48* ft; no roelmt d*Ut available in 1968. 
b ?reientea at "odd* altitudes tiouu original modal mi in os* unit*. 

°Preaoura above 426,504 ft la artrapoleted. 
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TABLE U. - HEEBCEP-aH L: SESAPOIAXSD A 3 MO 0 FEBBIC USHBL (REP. 3 ) 
[Cctmsrted to snglnaering units.] 


Altltnie 
(abort 
sea le-rel), 
ft 


0 

32 , BOS 
65 , BIB 
96,424 
131,232 
164,040 

196,848 

229,656 

282,464 

29 S ,272 

320,060 

360,888 

395,696 

426,504 

459,312 

492,120 

524,928 

557,736 

590,544 

623,352 

656,160 

721,776 

787,592 

853,008 

918,624 

984,240 

1 , 148,280 

1 , 312,320 

1 , 476,360 

1 , 640,400 


Physical properties 


ChaniQ&l ccBqooaifcLofT 


Tenp. far. 
raryiDg _/, 

or 


Pressure, Density, Ararage 7 aUm 
U>(P )/»1 ft lb(K)/cu ft aoltculsr percent 

wrtgfct, Og 


Ib/lb-aole 


2091 

577 

115 

26.4 

6.67 

1.87 

5 . 25 CXT 1 

1 . 21 XUT 1 

2 . 13 X 10" 2 

3 . 4 ®< 10" 3 

6 . 34 XUT* 

1 . 33 X 00 "* 

4 . 12 X 10- 5 

<=2.56X10-3 

9 . 93 X 10 "® 
6 . 3 ®ao-s 
4 . 64 XL 0" 6 
3 . 47 X 30 "® 
2 . 63 XL 0" 6 
2 . 05 X 10- 6 
1 . 62 X 10 “® 

1 . 06 X 10 "® 

7 . 21 X 10 "' 


2 . 44 X 10 "! 

l-UXlO " 7 

4 . 76 XUr® 

l.BTXUr® 

8 .BSX 10" 9 


7.463X10““ 

2 . 596 X 10 "® 

5 . 670 XLO " 3 

1 . 164 X 10" 3 

2 . 656 X 10 "* 

7.162XL0- 5 

2.2C6XL0" 5 

5,966X10"® 

1 . 224 X 10 - ® 

1.97000"* 

S^SSXlO ” 8 

5.696XHT 9 

1 . 337 XUT 9 

5.4S4X10" 10 

1.S0KX1O" 10 

7 .S 16 X 10" 11 

4 . 666 X 10 ”^- 

3.09000"^ 

2 . 11 S< 10 ' 1 - L 

1.574X10-11 

l.lBlXLO-H 

4 7.046XI0"?? 

4 . 576 X 30 "!® 

3 . 017 X 10 "!“ 

1.995X10” 12 

1 . 3 COO”!® 

s.saxur! 3 

2 . 53 XLD "” 3 

8 . 96 X 10 "!* 



Tbad.cal canqpon ltd. on. abcrvu - 426,604 ft 1 b 'batod on tboanrtlaal onion! gfrffgw 
^Presented at " 044 " altitudsa because or lgtml nodel me in age units. 
^Pressure aixTre 426,504 ft is eadj aj polat efl . 

^Density sborra 721,776 ft Is ortrspolatoi. 
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TABLE HI. - THE ARDC MODEL 


Altitude Physical properties 

(above 

sea level), Tamp, far Pressure, Density, 
ft varying JS, lb(F)/sq. ft lb(M)/cu ft 

Off 


246,075 

262,464 

295,272 

328,080 

360,888 

393,696 

426,504 

459,312 

492,120 

574,175 

656,160 

738,225 

820,250 

902,275 

984,240 

1,066,272 

1,148,280 

1,230,315 

1,312,320 

1,394,355 

1,476,360 

1.558.400 

1.640.400 


822.3 

922.4 

1024.4 

1127.5 

1231.0 

1334.7 

1438.1 

1541.3 

1544.1 

1746.5 

1848.4 

1949.7 

2050.4 


5.95X10”*; 

2.55X10"“ 

4.71XL0' 3 

9.66XL0” 4 

2.47X10" 4 

7.53XX0" 5 
c 2 .62X10"® 

l.nxio " 5 

5.63X10- 6 

1.58X10”° 

6.13X10” 7 
2.73X10” 7 
1. 34X10”! 
7.21XLO” 0 
4.11X10"® 

2.47X10”® 

1.55X10-S 

1.01X10”® 

6.81X10“° 

4.71XL0” 9 

3.53X10"® 

2.42X10”® 

1.78X10" 9 


5.095X10" 8 
1. 552X10”® 
2.494X10" 7 
4. 446X10”® 
9.920X10" 9 

2.670X10” 9 

B.109X10" 10 

2.682X10”!° 

l.rnxio - 10 

2.153X10” 11 

6.926X10" 12 

^.esoxio”!? 

1.177X10"" 

5.663X10”" 

2.936X10" 13 

1.622X10" 13 

9.39SX10" 14 

5.713X10"! 4 

3.582X10" 14 

2.328X10" 14 

1.556X10" 14 

1.071XL0”^- 4 

7.471XL0” 15 


a Estimated iron, graph in ref. 5. 

^Presented at "odd" altitudes to match other tables . 


^Pressure above 426,504 ft is extrapolated. 
^Density above 721,776 ft is extrapolated. 
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TABLE 17. - VERTICAL DISTRIBUTION OF ATOMIC OXYGEN AND OZONE 

(REF. 3) 


Altitude , 

Ozone 


Atomic oxygen 

ft 





Concentration 

Volume 

Concentration 

Volume 



of 0-, 

percent 

of 0, 

percent 


molecule/cc 


molecule/cc 


49,212 

1.0X10 11 

2. 44X10“ 5 





65,616 

2.0X10 12 

1.05X10 -4 



98,424 

6 . 7X10 12 

1.72X10" 3 

1.7X10 9 

4. 36X10" 7 

131,232 

4.8X10 11 

5.39X10” 4 

2.0X10 10 

2.25X10 -5 

164,040 

4.3X10 10 

1.79X10 -4 

4.4X10 10 

1.83X10" 4 

196,848 

7.0X10 9 

9.46X10" 5 

1.GX10 11 

1.35X10" 3 

229,656 

1.2X10 9 

6.00X10" 5 

1.9X10 11 

9.50X10" 5 

262,464 

7.6X10 7 

1.85X10' 5 

2.7X10 11 

6.59X10" 2 

295,272 

5.0X10 6 

7.58X10" 6 

3.3X10 11 - 

0.50 


TABLE V. - UPPER AND LOWER LIMITS OF ATOMIC OXYGEN RECOMBINATION 


IN NOZZLE OF RAMJET ORBITING AT 328,080 FEET IN ATMOSPHERE 

OF TABLE I 

[Nozzle length, 100 ft; frozen-compression diffuser.] 


Recombination- 
rate coefficient, 

kj, 

/ g-mole Y" 2 -1 

[msr) sec 

Near-equilibrium limit 

Near- frozen limit 

V 

g-mole 0 
liter 

Minimum 

mass 

percent 

recombined 

*0' 

g-mole .0 
liter 

Maximum 

mass 

percent 

recombined 

10 11 

5.3X10" 2 

0 

3.2X10" 6 

8 

10 12 

5. 3X10" 3 

0 

l.OXlO" 6 

32 

10 13 

5.3X10" 4 

0 

2. 7X10“ 7 

57 

10 14 

5.3X10 -5 

0 

6.6X10 -8 

78 

10 15 

3.7X10" 6 

4 

1.4X10“ 8 

91 

10 16 

1.2X10“ 6 

29 

2.7X10 -9 

97 
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Figure 2. - Variation of atmospheric energy available for pro- 
pulsion with altitude. - 


o elf 
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Figure 4. - Comparison of energy available with energy required 
for sustaining ionosphere orbits. Wall temperature, 1000° R. 
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Static temperature 



(a) Variation of tenperatura and preaaure. 

Figure S. - Thereod/namio ejelo for ionosphere raajet orbiting at 323,080 feet. 
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(b) Variation of enthalpy and preBsure . 

Figure 5. - Continued. Tharuodynasdo oyole far Ionosphere ramjet orbiting at 328,000 feet. 
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Figure 5. - Concluded. Thenaodynamic cycle for ionO' 
sphere ramjet orbiting at 328,080 feet. 
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Figure 7 . - Effect of ratio of inlet to exit area on ideal over-all engine efficiency 
for ramjet cycle of figure 5. 
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Figure 8. - Effect of ratio of inlet to exit area on thrust and drag for ramjet cycle 
of figure 5. Engine length, 100 feet; inlet radius, 20 feet. 




Net thrust, F, lb(F) 


NACA TN 4267 



Ratio of inlet to exit area, Aj/Aj 




4780 



Mass percent recombined 



--4 

H 


Altitude, ft 


Figure 11. - Variation of external equilibrium wall temperature with altitude. Engine length, 100 feet; nacelle 
half -angle, -8.6°; sunny aide. 




